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Abstract. For a class of one-dimensional mass transport models we present a simple 
and direct test on the chipping functions, which define the probabilities for mass to 
be transferred to neighbouring sites, to determine if the stationary distribution is 
factorized. In cases where the answer is affirmative, we provide an explicit method 
for constructing the single-site weight function. As an illustration of the power of this 
approach, previously known results on the Zero-range process and Asymmetric random 
average process are recovered in a few lines. We also construct new models, namely a 
generalized Zero-range process and a binomial chipping model, which have factorized 
steady states. 
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In a recent publication [I], we investigated a class of mass transport models on a 
ring (one-dimensional lattice with periodic boundary condition). The class encompasses 
both discrete and continuous masses and discrete and continuous time dynamical rules 
for the transfer of mass between neighbouring sites. The Zero-range process (ZRP) 
[3 Ej and Asymmetric random average process (ARAP) |3J E] correspond to special 
cases of this class. We derived a necessary and sufficient condition for the steady state 
to factorize which dictates an appealingly simple, yet general, form for the chipping 
functions (y?(/x|m) in equation ^ below) which define the probabilities with which mass 
is transferred from site to site. Though the form of this condition might be elegant, it 
is an "implicit test" for the chipping functions. 

In the present work we formulate a simple explicit test for the chipping functions, 
as well as a straightforward method for constructing the single-site weight function 
(products of which form the steady state distribution). In this sense, we provide a 
complete solution to this class of mass transport models. This note may be considered 
as a sequel to [I], where the reader will find the motivation for the model and details of 
the previous analysis. 

Following the notation of [J], our system consists of masses m« at site i = 1 ... L on 
a ring. At each time step, a mass /i; (drawn from a distribution <y2(/x,|mj)) is 'chipped off' 
mi and moved to site i + After long times, the system settles into a time- independent 
state, with distribution (unnormalized probabilities) F(mi, Wx). The central result 
of pQ (equation 15) is a necessary and sufficient condition on (p(jj,\m) for F to factorize: 
F = /(mi).../(mj, namely, if and only if ip(fi\m) can be expressed in the form 

[v * w\ (m) 

where v and w are two non-negative functions. Further, f(m), the single site weight is 
given by the convolution in the denominator, i.e., 

f (m) — [v * w] (m) = / d/x v (/x) w (m — /x) . (2) 

Jo 

Typically, a mass transport model is motivated by a specific chipping function ip. 
As a result, it is not particularly easy to see if it satisfies condition (P), the form of 
which, though elegant, is "implicit" . We now take the next simple step and turn this 
into an "explicit" test for <p. As in pQ, it is clearer to regard the two variables in (p as 
/x and a = m — /x. Note that the condition for factorization in Eq. ([TJ is equivalent to 
asking if 



d_ 
<9/x 



d_ 



ln^(/x|/x + a) (3) 

is a function of /x + a alone (i.e., m, and no other dependence /x, say). Thus, if we have 
factorization, the quantity in Eq. (j^J) is a function of m alone and let us call it h(m). 
That this is a necessary and sufficient condition equivalent to (0) can then be seen by 
integrating @ with respect to /x and a, which yields (p of the form (0). Moreover, 
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integrating twice with respect to m and then exponentiating, we have explicitly 



/ (m) = exp 



dm' / dm"h (m h 



(4) 



Note that there are two arbitrary integration constants in this construction, leading to 
an overall amplitude and exponential factor a m . These are precisely the two "degrees 
of freedom" we encountered while defining f(m) in pQ. 

This formalism thus provides a simple test @ to see if any specific mass transport 
model admits a factorized steady state and, if so, a recipe for the associated single-site 
weight (jlj). As an example of the applicability of (JHJ) and (jH), we consider the case of 
the ARAP [3J E] which is a class of mass transport models defined by a chipping 
function of the form (p(fj,\m) = tp(fj,/m)/m i.e. a random fraction r = fi/m chips off at 
each update. In this case ® gives 

.IJHS.^.J. (5) 

m z dr m z m z 

where 

^ ^ ip{r) dr ^ ^ 

Thus, the condition that Q depends on m alone implies h(m) oc 1/m 2 and 

— [r(l — r)K{r)} = constant (7) 
the solution of which is 

^(r) = Cr p {l -r) q , (8) 
where the constant C is fixed by the normalisation condition 

drip(r) = l. (9) 
This condition yields 



o 



r p(-\ — r \i 

M r ) = ^_ (10) 

VK ' B(p + l,q + l) K J 

where B(p+ 1, q + 1) is the usual Beta function and p,q > —1. The function h(m) given 
by is from © 

h(m) = \ \~ [r(l - r)K(r)] + l) = 1+p + q (n) 
[ dr J m z 

and, integrating twice (j3J), we retrieve the result: f(m) oc m 1+p+q . Thus fllOj) is the 
most general chipping function of r that gives rise to a factorized steady state. This 
proves in a direct way a result of jHj- 

For models with discrete masses such as the ZRP, the derivatives above become 
differences. To be careful, let us write the chipping rates as 

oo n 

ip(n\m) = 2j ^2 V^^G" - t)8{m - n) . 
n=l e=o 
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The factorization condition for continuous masses in Eq. 
analogue, 



<Pi,-. 



VlW r , 



fn 



has an equivalent discrete 



(12) 



where / n 's are single-site weights. The factorization test can now be rephrased in terms 
of the cross ratio 

.n 



R(£,n) = 



(13) 



^l+l,n+l^l,n+l 

defined when all of the (p's are positive. If the factorization condition in Eq. (|12j) 
holds then the cross ratio R is clearly independent of I and depends only on n. Thus, 
a factorized steady state exists if and only if R(l, n) is independent of I. In this case 
R{n) will be given in terms of the single-site weights f n in 



f(m) = ^2f n 5(m - n) 



n=l 



as 



R{n) 



f 2 

■J n 



+ 1 



fn fn+2 

which yields the recursion 

fn+2 _ 1 fn+1 

f n+1 R(n) f n 
Iterating (fTHjl yields 

fn+2 



1 

n ~ru 



m=0 



R(m) 



fa 



fn+1 

which one can iterate again to obtain 

"n-2 



(14) 



(15) 



(16) 



fn — {fo) 



n 

m=0 



1 



R(m 



,n—m—l 



for n > 2 . 



Of course, this is the discrete version of "integrating twice" . The two arbitrary constants 
(overall amplitude f and exponential amplitude a n = (/i//o) n ) are again explicitly 
displayed. 

As an illustration consider the ZRP where only unit masses can chip off with 
probability u(n) where n is the mass at the site In this case, we see that there 
is only one cross ratio, namely for / = 0, 

u (n + 2) [1 - u (n)] 



R(0,n) 



(19) 



u(n+l) [l-u(n + 1)] 
Since this R is automatically "independent of we immediately recover the conclusion: 
the ZRP admits a factorized steady state. Further, it is straightforward to retrieve the 
/ n 's from (fTHj) as 



fn 



fo 



u[n) 



fi<l) 
fo 



n 

m=l 



1 — u(m) 



mm) 



(20) 
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which recovers the result of [Q originally derived by a more complicated approach in 
jU]. Note that we have again displayed the two arbitrary constants with (...) brackets. 

As a simple example of a new model with a factorized steady state we define the 
binomial chipping model, in which the mass at each site is discrete, m = 0, 1, 2, . . .. The 
chipping is specified by the following kernel, 

<Pi,n = (fjp\l-p) n ~\ (21) 

where < p < 1 is a parameter and £ — 0, 1, . . . n. One can understand the model by 
interpreting n as the number of unit masses at a site, each of which move independently 
with probability p at each time-step. One finds 

77-1-2 

R(l,n) = ^TT ( 22 ) 
n + 1 

for all t. Since manifestly this is independent of t the steady state factorizes and (fTHj) 
yields 

f - A ($fh- (23) 

Note that, apart from the "irrelevant" factors, the single site weight here is extremely 
simple: 1/n!. 

As a more involved example of constructing a new model with a factorized steady 
state, we consider a generalized zero-range process where mass chunks of size one or two 
can chip off at each time step with probabilities ui{n) and uzin) respectively. In this 
case we have two cross ratios 

fl(0, B )= ; il ^"'7 ll :'" ,,WI ,o for n>0 (24) 
Mi (n + 1) [1 - «i (n + 1) - « 2 (n + 1)] 

M" + '>"■("> for „>!. (25 , 
u 2 [n + 1) «i (n + 1) 

If we demand that the stationary state factorizes, then we must have R (0, n) = R (1, n) 
which reduces to 

Mn + - ^(n) - u 2 (n)) =A ^ ^ 
tti(n + l)wi(n) 

where A is a positive constant independent of n. In terms of the ratio 



Pn = , 27 

uAn) 



becomes 



, 1 - u\(n) , 



«i(n) 

Though this condition is in the form of a nonlinear recursion, we can linearize it by 
changing variables to y n via 

Pn = A V -^ . (29) 

Un 
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Then (J28j) becomes a linear second-order recursion for the y's: 

1 — u\(n) 



Vn+l 



uiin 



-y n - Ay, 



n-l 



(30) 



with initial conditions yo = and y\ = 1. Thus we can find y n in terms of arbitrary 
Ui(n) and A. At the same time, these ?/'s also fix the "allowed" chipping rates for two 
mass units, U2{n), in terms of ui(m), m < n and one free parameter A, through the 
relations (|27j) and (|2?ty : 



Aui(n) 



2/n 



Of course, the single-site weights can also be found using (fTHj) and 



/n 



fo 



«i(n) - u 2 (^) 
which may also be written in terms of y n as 

7i«i(i)V y n 



/o 



n 

n 1 



m=l 



u\{m) — U2{m) 
Ui(m) 



fn — fo 



fo 



ui(n) 



(31) 



(32) 



(33) 



Finally, we remark that the constructive method presented here easily adapts to the 
case of continuous time (or random sequential) dynamics specified by the rates per unit 
time, 7(/i|m), at which mass [i chips off mass m. It is shown in that the necessary 
and sufficient condition for a factorized steady state becomes 

x(fi)w(m — ji) 



7(^|m) 



w(m) 



(34) 



in which case the single-site weights become f(m) = w(m). Thus the test for and 
construction of the factorized steady states is identical to that of the discrete time case 
with ip(fj,\m) replaced by 7(/i|m). 



Acknowledgments 

This research is supported in part by the US National Science Foundation through 
DMR-0088451 and DMR-0414122. 

[1] M.R. Evans, S.N. Majumdar, and R.K.P. Zia, J. Phys. A 37, L275 (2004) 

[2] F. Spitzer Advances in Math. 5 246 (1970) 

[3] M.R. Evans, Braz. J. Phys. 30, 42 (2000) 

[4] J. Krug and J. Garcia, J. Stat. Phys. 99, 31 (2000). 

[5] R. Rajesh and S.N. Majumdar, J. Stat. Phys. 99, 943 (2000). 

[6] S.N. Coppersmith, C.-h. Liu, S. Majumdar, O. Narayan, T.A. Witten Phys. Rev. E., 53, 4673 
(1996) 

[7] F. Zielen and A. Schadschneider, J. Stat. Phys 106, 173 (2002) 
[8] F. Zielen and A. Schadschneider, J. Phys. A 36, 3709 (2003) 
[9] M. R. Evans, J. Phys. A: Math. Gen. 30 5669 (1997) 



